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Abstract. This work reports a theoretical and experimental study of 4H-SiC JFET threshold voltage 
as a function of substrate body bias, device position on the wafer, and temperature from 25 °C 
(298K) to 500 °C (773K). Based on these results, an alternative approach to SPICE circuit 
simulation of body effect for SiC JFETs is proposed. 

Introduction 

Over the past few years, significant progress has been made in the development of increasingly 
capable very high-temperature integrated circuits (ICs) that could benefit a variety of aerospace, 
automotive, and energy production systems. For example, elsewhere at this conference it is reported 
that a 4H-SiC Junction Field Effect Transistor (JFET) IC with two levels of metal interconnection 
has demonstrated 500 °C electrical operation in excess of 3000 hours [1]. Increasingly accurate and 
accessible IC design and simulation tools for such extreme temperatures are clearly relevant to 
further technology development and application adoption. Towards this end, this work reports a 
theoretical and experimental comparison study of 4H-SiC JFET threshold voltage (Vt) as a function 
of substrate body bias (Vs) and temperature (T) from 25 °C (298K) to 500 °C (773K). 


Experiment 

Fig. 1 schematically depicts the cross-section of 
the long-channel (gate-length L (l = 6pm) 4H-SiC 
epilayer JFET structure studied. Key structural 
terms illustrated in Fig. 1 include gate and 
substrate pn junctions, n-channel epilayer 
thickness t c and doping Ndc, and sub-channel 
epilayer doping Nas- Experimental procedures 
used for device fabrication, pre-dicing wafer 
probe mapping of JFET V T , and high temperature 
packaging and testing are described in [1,2]. 

Theory. Under one-dimensional depletion 
approximation [3], the depletion depth x n on the 
n-side of a pn step-type junction with built-in 
voltage Vbi under applied pn junction bias Va is: 



Fig. 1. JFET Cross-sectional diagram. 
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Na and Nd are junction p-type and n-type layer doping densities, e s is 4H-SiC dielectric constant, q 
is the electron charge constant, and junction n-side depletion constant K n defined as: 



K = |f 2e sV n a \ 
n aJV q / VN d (N a +N d )/ 


( 2 ) 


When the JFET gate bias Vg is at threshold voltage Vt, the gate and substrate pn junction depletion 
depths extending into the n-channel (x n o and x n s) merge together to effectively remove carriers 
(electrons) from the entire doped n-channel thickness beneath the p + gate. At threshold therefore, 
the sum of the gate and substrate junction n-side depletions approximates the n-channel thickness t c : 

tc = XnG + Xns = K nG V(V biG - V T ) + K nS V(V bis - V s ) (3) 

where Vs is the applied substrate junction bias. Solving Eq. 3 for V T : 

Vt = V blG - (£)’ + (^) - (f^) 2 (V bls - V s ) (4) 

By defining Vto = Vt at Vs = 0, Eq. 4 can be re-expressed as: 

Vt = V T0 + (2g=s) (VVbis - V s - V S (5) 

Measured Results. Consistent with results 
from our previous study of a different 4H-SiC 
JFET IC wafer [2], the magnitude of Vt on this 
wafer increased proportional to the square of 
radial distance ( r ) from the wafer center (Fig. 2, 
blue). Secondary Ion Mass Spectroscopy (SIMS) 
analysis was performed (by Evans Analytical 
Group, www.eag.com) on two diced 3mm x 3mm 
die taken from differing radial wafer locations 
with significantly different measured JFET V T . 

The Fig. 3 SIMS results indicate that n-epilayer 
thickness variation is primarily responsible for 
the wafer position variation of Vt. Given 
position-invariant 1.0 x 10 cm" n-channel 
epilayer doping Ndc evidenced by Fig. 3 SIMS 
data, Eq. 3 can extract n-epilayer thickness t c as a 
function of r (Fig. 2, red), which like Vt exhibits 
y = a *r 2 + b general behavior. However, it should be noted that variation in p-type sub-channel 
doping N A s at/below SIMS aluminum detection floor (~ 3 x 10 cm" ) could slightly affect V T 
(~2 V, as seen in Fig. 4). A JFET diced from a different wafer location was custom-packaged for 
high temperature oven electrical testing [2], 

Comparison with Theory. Fig. 4 plots theoretical Vt vs. Vs curves (lines) calculated via Eq. 4 
using both Fig. 3 SIMS profiles assuming three different sub-channel doping densities (Nas =1,2, 
and 3 x 10 cm" ). Good agreement with V T experimentally measured on corresponding JFETs 
(plotted as symbols) is obtained. Fig. 5 compares measured (symbols) and theoretical (lines) V T at 
Vs = OV, -15V, and -25V from 25 °C to 500 °C from the high-temperature packaged JFET. The V T 
temperature dependence arises from junction built-in voltage Vbi G and Vi„s temperature dependence. 
Fig. 5 theoretical Vt approximations fall within 15% of measured V T , which is reasonable 
agreement considering uncertainty in Nas and t c for this device that was not SIMS profiled. 

SPICE Model. As SPICE [4,5] is the most broadly available electrical circuit simulation/design 
program, it is important to implement reasonably accurate and accessible SPICE models for 4H-SiC 
JFET electrical behavior. To ensure that such models can run on all variants of SPICE software, 
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Fig. 2. Measured Vt (blue) and extracted t c 
(red) vs. r from 25 °C wafer probe map data. 
The wafer center is Die 13,13. 
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Fig. 3. SIMS of JFET epilayers at two locations 
on the wafer. 



Fig. 4. Theoretical (lines, calculated via 
Eq. 4) and measured (symbols) JFET V T 
vs. Vs for die characterized by SIMS. 


this modeling should be implemented using baseline-version SPICE device models. However, the 
baseline SPICE JFET model makes no provision whatsoever for body bias effect. Therefore, the 
baseline SPICE NMOS model, which importantly includes body bias effect, is instead chosen for 
first-order SPICE modeling of 4H-SiC JFET electrical behavior. It is important to note that this 
NMOS modeling approach is only valid so long as gate and substrate junctions do not become 
forward biased (which would draw substantial gate/substrate junction current in a physical JFET 
not accounted for in the NMOS SPICE model). While this general approach was previously used 
for 6H-SiC JFET SPICE modeling [6], this previous work did not study body bias effect. Eq. 6 
shows the baseline SPICE NMOS body effect model formula for V T as a function of Vs in terms of 
baseline SPICE NMOS model parameters PHIB, VTO, and GAMMA [4,5], 


V T = VTO + GAMMA(Jl ■ PHIB - V s - V2 ■ PHIB ) 


( 6 ) 


Comparison of Eq. 5 to Eq. 6 readily reveals that 
they have the same mathematical form except for 
the linear Vs term residing at the end of Eq. 5. 
For SiC epitaxial JFET IC designs of interest (i.e., 
those that approximate Fig. 1 with |Vs| < 50V), 
the linear Vs term is numerically the smallest 
term of Eq. 5. 

Table I shows two proposed mappings for 
substituting 4H-SiC JFET physical parameters 
directly into the Eq. 6 SPICE NMOS body effect 
model. Parameter Map #1 is obtained by directly 
mapping Eq. 5 terms into obviously 
corresponding Eq. 6 terms while completely 
neglecting the linear Vs term at the end of Eq. 5. 
Map #2 features a mathematically adjusted 
GAMMA that approximates the impact of the 
linear Vs term at the end of Eq. 5 within the form 
of Eq. 6. The Vt vs. Vs plots of Fig. 6 compare 
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Fig. 5. Theoretical (lines) and measured 
(symbols) temperature dependence of 
JFET V T at V s = 0V, -15V, and -25V. 


Table I. Parameters for SiC JFET body effect modeling. 
SPICE NMOS JFET Parameter JFET Parameter 


Parameter [4] 

Map #1 


Map #2 

VT0 

Vto 


Vto 

PHIB 

0.5 V blS 


0.5 V blS 

GAMMA 

2t c K n s 

K nS 

K^g 

(2t c - K nG ) 



these two SPICE Eq. 6 approximations against the full 
Eq. 5 body effect theory for the two Fig. 3 SIMS 
devices. As seen in Fig. 6, improved agreement with 
Eq. 5 is obtained for Parameter Map #2. 

It is important to note that the NMOS model will 
crash baseline-version SPICE software when the SPICE 
TEMP (temperature) parameter exceeds -300 °C. This 
fact precludes direct use of TEMP=500 input with 
baseline-version SPICE to simulate 500 °C JFET 
behavior starting from a 27 °C (SPICE default TEMP) 

NMOS model [6]. Therefore, a different/unique “500 °C” NMOS model (with parameters adjusted 
to match 500 °C JFET behavior) must be run with TEMP=27 to carry out baseline -version SPICE 
simulations of 500 °C circuits. This approach to SPICE modeling the full characteristics of 4H-SiC 
JFETs and circuits for temperatures up to 500 °C will be presented in more detail in future work. 


Fig. 6. Comparison of Vt vs. Vs 
calculated for Fig. 3 SIMS JFETs at 
T= 25 °C. Black solid lines are Eq. 5, 
while colored/dashed use the Eq. 6 
with parameter maps listed in Table I. 


Summary 

This work has derived and experimentally verified equations for modeling the body/substrate bias 
effect of 4H-SiC n-channel epilayer JFETs over a wide range of temperature and substrate bias. 
Parameters for using the baseline-version SPICE NMOS transistor model for simulating 4H-SiC 
JFET body bias behavior have been described. This study has also ascertained that n-channel 
epilayer thickness variation (not doping variation) is responsible for the observed increase in V T as a 
function of the radial distance from the center of the experimental 76mm diameter wafer. 
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